c->cd(4);
h4->ProfileX("_pfx1", 0, 400);
_pfx1->Draw("same");
TF1 *rt = new TF1("rt", "pol3");
rt->SetParLimits(0, -0.2, 0.9);
_pfx1->Fit(rt, "QBR", "", -20., 350.);
for(int i = 0; i <= 3; i++){ cout << rt->GetParameter(i) << (i<3?", ":"\n"); }
cout << "r(0) = " << rt->Eval(0.) << " r(162) = " << rt->Eval(162.) << endl;

c->cd(5);
h5->ProfileX("_pfx2", 0, 400);
_pfx2->Draw("same");
TF1 *rt2 = new TF1("rt2", "pol3");
rt2->SetParLimits(0, -0.2, 0.9);
_pfx2->Fit(rt2, "QBR", "", -20., 350.);
for(int i = 0; i <= 3; i++){ cout << rt2->GetParameter(i) << (i<3?", ":"\n"); }
cout << "r(0) = " << rt->Eval(0.) << " r(162) = " << rt->Eval(162.) << endl;

c->cd(6);
h6->ProfileX("_pfx3", 0, 400);
_pfx3->Draw("same");
TF1 *rt3 = new TF1("rt3", "pol3");
rt3->SetParLimits(0, -0.2, 0.9);
_pfx3->Fit(rt3, "QB", "", -20., 350.);
for(int i = 0; i <= 3; i++){ cout << rt3->GetParameter(i) << (i<3?", ":"\n"); }
cout << "r(0) = " << rt->Eval(0.) << " r(162) = " << rt->Eval(162.) << endl;

